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Abstract. 

We study the nonequihbrium properties of directed Ising models with non conserved 
dynamics, in which each spin is influenced by only a subset of its nearest neighbours. 
We treat the following models: (i) the one-dimensional chain; (ii) the two-dimensional 
square lattice; (iii) the two-dimensional triangular lattice; (iv) the three-dimensional 
cubic lattice. We raise and answer the question: (a) Under what conditions is the 
stationary state described by the equilibrium Boltzmann-Gibbs distribution? We show 
that for models (i), (ii), and (iii), in which each spin "sees" only half of its neighbours, 
there is a unique set of transition rates, namely with exponential dependence in the 
local field, for which this is the case. For model (iv), we find that any rates satisfying 
the constraints required for the stationary measure to be Gibbsian should satisfy 
detailed balance, ruling out the possibility of directed dynamics. We finally show 
that directed models on lattices of coordination number z > 8 with exponential rates 
cannot accommodate a Gibbsian stationary state. We conjecture that this property 
extends to any form of the rates. We are thus led to the conclusion that directed 
models with Gibbsian stationary states only exist in dimension one and two. We then 
raise the question: (b) Do directed Ising models, augmented by Glauber dynamics, 
exhibit a phase transition to a ferromagnetic state? For the models considered above, 
the answers are open problems, to the exception of the simple cases (i) and (ii). For 
Cayley trees, where each spin sees only the spins further from the root, we show that 
there is a phase transition provided the branching ratio, q, satisfies > 3. 
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1. Introduction 

A long-standing theme of the statistical physics of nonequilibrium systems is the 
question of the nature of their stationary state, and in particular of the existence of 
phase transitions at stationarity. By nonequilibrium systems, we mean systems whose 
dynamics is non reversible, i.e., such that the dynamical rules defining the model do 
not obey detailed balance. For example, systems which are driven, e.g. submitted to a 
field, do not obey detailed balance. The nature of the stationary state of driven systems 
has been well investigated for a number of models, among which are the Zero Range 
Process [1], the Asymmetric Simple Exclusion Process [2], or the KLS model [3, 4]. 

In contrast, the same questions have not yet been addressed in a comprehensive 
way for directed Ising models with non conserved dynamics. In one dimension, for 
example, a spin is influenced only by its right neighbour, while for a two-dimensional 
square lattice, for example, each spin is influenced only by the neighbouring spins above 
("North") and to the right ("East"), and not by the spins below and to the left. The 
asymmetry of the interactions means that the detailed balance condition is not obeyed, 
and therefore there is no conventional equilibrium. Again one speaks instead in terms of 
a stationary state. However, the latter differs from the stationary state that prevails for 
driven systems in that it docs not carry a macroscopic current. In the present work we 
endeavour a study of these questions, namely the nature of the stationary state measure 
and the existence of phase transitions, for directed Ising models with single spin-flip 
dynamics. 

Usually, when the dynamics of a model does not fulflU detailed balance, one expects 
its stationary measure to be non-Gibbsian. Consider for instance the Ising chain with 
conserved dynamics, and totally asymmetric rules: only the H — bond is updated, not 
the — h bond. As pointed by Katz, Lcbowitz and Spohn [3], while generically the 
stationary measure of this model is non-Gibbsian, for special choices of the rates defining 
the dynamics, the stationary measure is the same as that of the model with symmetric 
dynamics, i.e., the equilibrium Boltzmann-Gibbs distribution. The one-dimensional 
KLS model, which can also be viewed as the hopping of one species of particles with 
exclusion on a lattice, has been further analysed and generalised to the case of two 
species of particles hopping with exclusion [4], where, again, the stationary measure 
is the same whether the system is driven or not. To date, for spin models with non- 
conserved dynamics, no systematic analysis of this kind has been performed. The main 
purpose of this paper is to provide such an analysis, for the Ising model, in one, two, 
and three dimensions. 

We treat the following models: (i) the one-dimensional chain; (ii) the two- 
dimensional square lattice; (iii) the two-dimensional triangular lattice; (iv) the three- 
dimensional cubic lattice. The question to be answered is how to choose the transition 
rates, for the single spin-flip dynamics, in such a way that they lead to a stationary state 
with the equilibrium measure (whether or not the dynamics is directed)? The method 
consists in writing the master equation at stationarity, then imposing the equilibrium 
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measure [4]. One thus finds a set of constraint equations on tlie rates. Additional 
constraints come from the choice of dynamics (with or without spin symmetry, directed 
or not). For the undirected models (i.e., with symmetric dynamics) we find, as expected, 
that any rates leading to the Gibbs measure satisfy detailed balance. For the directed 
case, we show that for models (i), (ii), and (iii), in which each spin "sees" only half 
of its neighbours, there is a unique set of transition rates, namely with exponential 
dependence in the local field, for which this is the case. By local field we mean the 
restricted field felt by the fiipping spin. For a variant of model (ii), where each spin sees 
only its West, North and East neighbours, the property still holds but the set of rates 
have no longer a simple dependence in the local field. For model (iv), we find that any 
rates satisfying the constraints required for the stationary measure to be Gibbsian should 
satisfy detailed balance, ruling out the possibility of directed dynamics. We finally show 
that directed models on lattices of coordination number z > 8 with exponential rates 
cannot accommodate a Gibbsian stationary state. We conjecture that this property 
extends to any form of the rates. We are thus led to the conclusion that directed models 
with Gibbsian stationary states only exist in dimension one and two. 

The second purpose of this paper is to investigate whether directed Ising models, 
augmented by Glauber dynamics, can exhibit a phase transition to a ferromagnetic state. 
This question was motivated by a recent paper [5], where Lima and Stauffer considered 
Ising models with directed interactions in two to five space dimensions, evolving via 
Glauber dynamics [6]. They find in particular that for the two-dimensional square 
lattice where each spin is infiuenced only by the North and East neighbouring spins 
there is no ferromagnetic phase transition, in contrast to the usual equilibrium Ising 
model which has a second-order phase transition in two dimensions. In the present 
work we show how the absence of a phase transition arises both for the directed Ising 
chain and for the two-dimensional directed square lattice considered in [5]. We then 
investigate Cayley trees, with the influence directed from the tips towards the root, and 
show that there is a phase transition for coordination number z > A but not for z — ?>, 
i.e., for branching ratio > 3, since q = z — 1. 

The outline of the paper is as follows. We begin, in the following section, by 
giving an extensive analysis of the one-dimensional case. We show in particular that, 
for the directed Ising chain, there is a unique set of rates which produce stationary 
states described by the Gibbs measure. We then address in section 3 the case of two- 
dimensional and three-dimensional Ising models on regular lattices. Section 4 is devoted 
to another approach, which allows to determine which lattices, with given coordination 
number, can sustain directed models with rates of exponential form. We then proceed, 
in section 5, by studying directed models with Glauber dynamics. We show that for the 
two-dimensional square lattice there is no phase transition to a ferromagnetic state, the 
latter result being in agreement with the numerical data of [5]. The same result holds 
for the directed Ising chain. We show that for the Cayley tree model, a phase transition 
occurs if the branching ratio is greater than or equal to three. Some generalizations 
are given in the Appendix. Section 6 concludes with a discussion and summary of the 
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results. 



2. Equilibrium Measures for Nonequilibrium Stationary States 

In this section we pose the question: "Under what conditions are the stationary states 
of Ising models described by the Boltzmann-Gibbs distribution?" . We start by the one- 
dimensional Ising model. We consider higher dimensional models in sections 3 and 4. 

2.1. The method illustrated on the case of the Ising chain 

We derive a set of constraints between the transition rates that need to be satisfied in 
order for the stationary state measure to be Gibbsian. For symmetric dynamics, the 
constraints are satisfied if and only if the rates satisfy detailed balance. For directed 
dynamics, however, the constraints uniquely determine the rates (up to an overall 
timescale) . 

The energy of configuration C = {(7i, . . . , cTj, . . . , un} is given by 

E{C)^-JY.<Tm+i, (2.1) 

i 

and we choose periodic boundary conditions. The dynamics consist in flipping a spin, 
chosen at random, say spin i, with a rate w{Ci\C) corresponding to the transition between 
configurations C and Ci — {(Ti, . . . , — (jj, . . . , cr^v}- The change in energy due to the fiip 
reads 

/^E = E{Ci) - E{C) = 2a,h, = 2a, J{a,_^ + C7i+i), (2.2) 

where we denote the local field acting on spin i by hi. At stationarity, the master 
equation expresses that losses are equal to gains, and reads 

P{C) Y.<C,\C) = Y.w{C\Ci)P{C,), (2.3) 

i i 

where, by hypothesis, 

P{C) oc e-^-^(^). (2.4) 

After division of both sides by the weight P{C), eq. (2.3) can be rewritten as 

Y^widlC) - wiC\Ci)e-^^^ = 0, (2.5) 

i 

which is also, in the present context, the general form of the stationary master equation 
in any dimension. 

There are eight (2^) possible rates w{Ci\C) denoted by w((Ti_i CTj (Tj+i), or for short 
by 'Ji'c7i_i CTi+i if = +1) or by Wo-j-i o-i+i if — ~1 (see Table), corresponding to the 
fiipping of spin (jj for the eight possible motifs (in the present case, the triplets appearing 
in the Table): 

(7j_i CTj (Tj+i ^ (7j_i (-CTj) (7j+i. (2.6) 
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Table 1. List of moves and corresponding energy difference AiJ (eq. (2.2)) for the 
Ising chain. 



Move 



+ + + 
+ + - 

- + + 

- + - 



+ 
+ 



+ 



+ 



Rate 



w 



AJ 



-AJ 



Move 



+ 
+ 



+ 



+ 



+ + + 
+ + - 

- + + 

- + - 



Rate 



w 



-AJ 



AJ 



Let A^o-j_i iTi o-j+i be the number of such motifs in configuration C, now considered as 
fixed. The balance equation (2.5) is finally recast as 

N+++{w++ - e"^^ W++) + iV++_(w+_ - w+J) + iV_++(w_+ - 
+A^_+_(w__ - e^-^ w__) + N+_+{w++ - w++) + iV+__(w+_ 



+A^__+(i(;_+ - W-+)+N {w-- - e-^^w^J) = 0, (2.7) 

where we denoted the reduced coupling constant (3 J by K. A shorter notation for the 
same equation is 



^iV(/.„) [«;(/.„) -«;(/.«)e-^(^(^«)-^('^")) 



a=l 



(2.^ 



where N{^a) is the number of occurrences of the motif //^ in configuration C, and Jia is 
the motif obtained after flipping the central spin in the motif. 

Equation (2.8) should be satisfied for any given configuration C. Since the quantities 
inside the brackets are couples of rates related by detailed balance conditions, the 
left side of (2.8) vanishes identically as soon as these conditions arc imposed on the 
rates. However imposing detailed balance is not a priori necessarily the only way of 
solving (2.8), because the N{na) are not independent quantities. One should therefore 
express the latter on a basis of independent quantities, in the present case correlators 
defined as follows. Let us denote by = (l+(Ti)/2 the indicator variable for the presence 
of a + spin on site i. Then, for example. 



N 



N 



N 



8 

where, for example. 



i=l ^ i=l 

1 + Ci + C2 — C3 + Ci2 — Ci3 — C23 — C123), 



(2.9) 



N 



N N 

— ^CTi+i, C12 = — X^c^iCTm, 

1=1 i=l 



ci^ — }^<^i, C2 = —}^(Ji+i, C12 = — } ^(Jiai+i, (2.10) 
and so on. Translational invariance imposes the following identities between correlators, 

Cl = C2 = C3, C23 = C12. (2.11) 
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Rewriting the master equation (2.8) in terms of the remaining independent correlators 
yields 

ei + Ci 62 + Ci2 63 + Ci3 64 + C123 65 = 0, (2-12) 

where the coefficients 61, . . . , 65, are linear combinations of the expressions appearing 
in the brackets of eq. (2.8). This equation should be satisfied for any configuration 
C, hence the coefficients of the correlators must vanish identically. This provides five 
constraint equations, Ci = . . . = 65 = 0, not all independent, on the eight rates w^cr ^-nd 
w^a- Reducing this set of equations yields the final constraint equations: 

w++^Q^^w++ (2.13) 
= e-^-^w__ (2.14) 
+ lO-^ = + (2-15) 

The first two equations are detailed balance conditions. They involve rates for moves 
corresponding to a non-zero value of the local field J(crj_i + (Jj+i) = ±2J. The third 
equation involves the rates for moves that do not imply a change in the energy, i.e. with 
zero local field (motion of a domain wall). 

The eight rates fulfilling the constraints (2.13), (2.14, (2.15) therefore depend on 
five independent parameters. A general expression of these rates is 

w{(Ti 0-2(73) = ^|l + A(7i - A'(72 + ^2^ - (73 + 5 (7i(73 + e(7i(72 

-(7(1 + (5) + e) (72(73 - A'(7i(72(73|, (2.16) 

where A, A', o;, 5, e are the five independent parameters, while 7 is given by: 
e^^-1 
e^^ + l 

and (7i, (72, (73 are respectively the left, central, and right spins. These rates correspond 
to the generic situation where the central spin is influenced unequally by the left and 
right spins. 

Hereafter we restrict the discussion to the case of rates invariant by spin symmetry, 
that is 

w^^ = ty-o-(T- (2.18) 

In this case, the constraints (2.13), (2.14, (2.15) reduce to a single equation, 

w--^e^^w++, (2.19) 

hence three rates remain unknown, namely Setting A = A' = 0, in order 

to satisfy (2.18) we obtain 

w{ai (T2 (T3) = - (1 + 5 (Ti(T3 + e ai(72 - (7(1 + 5) + e) a2a3) , (2.20) 

depending now on the three independent parameters, a, 5, e. 

We now address two extreme cases of interest, where the dynamics is symmetric, 
then when it is totally asymmetric, or directed. 



7=;4^-r=tanh2X, (2.17) 
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2.2. Symmetric dynamics 

The dynamics is spatially symmetric, or undirected, if the left and right spins have an 
equal influence on the central spin. This right-left symmetry imposes therefore that 

w+_ = (2.21) 

This fixes one of the remaining unknown rates. By (2.21) we obtain 2e + 7(1 + 5) = 0, 
and therefore 

w{a^ (72 (73) = I (^1 + (5 (71(73 - ^7(1 + ^) <72((7i + (Js)) , (2.22) 

which is the most general form of the rates for symmetric dynamics, originally proposed 
by Glauber (see Appendix of [6]), which now depend on two parameters. This form 
satisfies the detailed balance conditions (2.19) and (2.21). 

Simplified forms of the rates are obtained by fixing the parameters a and 5. 
Examples of possible choices are the Glauber (or heat bath) rule: 

w{a, (72 (73) = ^ (1 - (72 tanh(X((7i + (73))) = ^^A^qry, (2.23) 
or the Metropohs rule: 

w(ai (72 (73) = min (l, e"^^^) , (2.24) 
or yet the following rule: 

w{a^ (72 (73) = e-P^^/'' = e-Ka^^'^^+a^) _ (2.25) 

where AE', eq. (2.2), is the change in energy of the system due to the flip. 

2.3. Directed dynamics 

Let us now consider the case where the dynamics is totally directed, i.e., totally 
asymmetric. Assume for instance that spin Ui only looks to the right, hence that the 
rates only depend on the right neighbour of the flipping spin: 

W++ = w_+, w+_ = w (2.26) 

Carrying these conditions into (2.20), we obtain 5 = e = 0, and therefore 

a 

w{(Ti (72 (73) = - (1 - 7 (72(73) . (2.27) 

This expression can be equivalently obtained by suppressing all terms where the spin 
(7i appears in the general expression of the rates (2.16). We thus find a unique solution 
to the problem posed, up to the global time scale a. Fixing this scale by the choice 
a — 2 cosh 2K, we obtain the compact form 

«;((7if72(73) =e-'^'^^"% (2.28) 

while, with the choice o; = 1, we obtain the alternate form 

w((7i (T2 (73) = ^ (1 - (72 tanh(2i^a3)) . (2.29) 
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The directed dynamics is more constrained than the undirected one, as previously 
examphfied for conserved dynamics of two or three species of particles in [4]. No 
kinetically constraint models (i.e., with vanishing rates) can be devised in this case, 
in contrast with the case of the undirected dynamics which allows more freedom. 

Let us finally mention that, as a special case of the generic expression (2.20), where 
the central spin is unequally influenced by its neighbours, the following form of the 
rates has the virtue of interpolating between the undirected case (2.25), and directed 
case (2.28): 

w(ai as as) = e-'^Ka2{xa^+{^-x)as) ^ (2.30) 
with < a; < 1. 



3. Ising models on regular lattices 

In this section we consider two- and three-dimensional Ising models on regular lattices. 
The energy of configuration C = {a"i, . . . , a"j, . . . , ajv} is now given by 

E{C)^-J^aiaj, (3.1) 

{id) 

where (i,j) are nearest neighbours. We follow step by step the method used in the 
previous section for the determination of the transition rates leading to a stationary 
measure given by the Boltzmann-Gibbs distribution. We start by the two-dimensional 
Ising model on the square lattice. We proceed with the case of the two-dimensional 
triangular lattice, then of the three-dimensional cubic lattice. 



3.1. Constraint equations on the rates for the square lattice 

On the square lattice each spin has four neighbours. The change in energy when flipping 
the central spin, denoted by a, reads 

AE = 2aJ{aE + (TN + (Tw + (Js), (3.2) 

where a^ is the East spin, the North spin, etc.J The stationary master equation can 
be compactly written as in (2.8) with rates corresponding to the 32 motifs aascrNcrwC^s- 
We use the following abridged notations: 

= W5, ty++-+- = We, + = W7, w++ = ws 

= Wg, = Wio, «^+-+-+ = Wu, «^+-+— = Wi2 
W+ ++ = Wis, W+ +_ = Wi4, W+ + = ^15, W+ = ^16, 



I The notation ajv, where N stands for North, should not be confused with the notation for the spin 
with index N, size of the system. 
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and similarly, 

= Wl, W-+++- = W2, = W3, = tV4 

W-+-++ = W5, = We, W-+ — + = wr, w_+ = 

W +++ = Wg, W ++_ = Wio, W +_+ = Wii, W + = W12 

W ++ = W12, W +_ = Wu, W + = Wi5, W = WiQ. 

The method then proceeds as in the one-dimensional case. The identities imposed 
by translational invariance on the correlators are 

ci = C2 = C3 = C4 = C5, (3.3) 

C12 — Ci4, Ci3 = Ci5, C25 = C34, C23 = C45, (3-4) 

where the indices 1, 2, . . . , 5, correspond respectively to the central spin, the East, North, 
West, South spins. After reduction, the equations on the rates yield fourteen constraint 
equations. If one restricts the study to situations where the rates are invariant by spin 
symmetry, then the following additional constraints must be taken into account: 

Wl = Wie, W-2 = Wi5, W2, = ^14, W4 = Wi3, 

= W12, Wq = u^ii, W'j = Wio, ^8 = twg, 

Wg = ^8, Wio = W'J, Wii = Wq, Wi2 = Ws, 

1^13 = Wa, ^14 = Wz, Wl5 = W2, Wm = ^^^1, (3.5) 

which in turn reduce the former system of fourteen constraint equations to a system of 
six equations, as follows: 



Wq-Wq = 0, 
AK 

3 W2-W2- 

AK - /„,. AK - 



e^^ W2-W2 + e^^ W5- W5 = 0, 



e W3 ~ ws - [ws - c ws) = 0, 

iTe' 



e^^ W2-W2- (e^^ W3 - W3) + , , ^4k ('^'^ ~ "^7) = 0, 



AK - , 4K - -~ ( - \ r\ 

e W2-W2 + e W3-W3-— TTr(W4 - W4) = 0. 

The first two equations involve couples of rates related by detailed balance conditions, 
the following ones linear combinations of such couples. The rates now depend on ten 
independent parameters. 

3.2. Symmetric dynamics on the square lattice 

The dynamics is spatially symmetric, or undirected, if the left and right spins, or up 
and down spins have an equal influence on the central spin. For simplicity, we restrict 
the study to the most symmetric dynamics where the rates only depend on the number 
of up and down neighbouring spins of the central spin, i.e., rates which only depend on 
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the value of the sum cr^ + ctjv + cfw + ^s- 

^4 = We = ^7 = WiQ = Wii = = z, 

Ws = Wi2 = Wu — Wi5 = t. (3.6) 

Carried in the constraint equations above, we find 

t = e^^ y, wie = e^^'wi, (3.7) 

and z remains arbitrary. After fixing the overall scale of time the transition rates still 
depend on two independent parameters. Fixing these parameters leads to the following 
simplified forms of the rates, which are simple generalisations of the ID expressions, 

w{aaEaN(Twcrs) ata,nh{K{aE + (Tn + (Jw + crs))) , (3.8) 



wlcraEO-NO-wf^s) = min(l,e ^^^j , (3.9) 
w{aaEaNaw(Js) = e'^^^^/^ _ ^-Kaia^+aN+Tw+Ts) _ (310) 

3.3. Directed dynamics on the square lattice 

Consider first the case where the central spin is only infiuenced by the North and East 
spins. This means that 

Wi — W2 — W3 — W4 — X, 

W5 ^ We ^ wt ^ Ws ^ z, 

Wq = Wio = Wii = W12 = z', 

Wl3 = Wu = ^"15 = ^"16 = y- (3.11) 

Carrying these relationships in the constraint equations (3.6) leads to 

y = e^^ X, z' = z, z = e^^ x. (3.12) 

Out of the four unknown x, y, z, z', only one remains undetermined, i.e., there only 
remains one free parameter in the expression of the rates. The general expression of the 
rates satisfying these constraints reads 

w{a ue (Tn ctw crs) = f + ")'^ ^^^-^ ~ lcr{crE + (Tjv)) , (3.13) 

which is the unique solution of the question posed, up to the global timescale a. Fixing 
this timescale by the choice a — 2cosh^2i^, allows to recast (3.13) into the compact 
form: 

w{a ge (Tn (Jw (^s) = e-'-^'^('^^+^~). (3.14) 

Such rates violate detailed balance, but yet, as for the one-dimensional case, lead to a 
stationary state with Gibbs measure. This form of the rates for the two-dimensional 
directed Ising model on the square lattice where each spin sees only its North and East 
neighbours was first proposed in [7] (up to a seemingly missing factor 2 in the exponent) . 
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We also considered the more general case where the central spin is influenced by 
three of its neighbours, say, West, North, East. Introducing the simplifying notations 
wi, . . . , ^8, we have 

= W2 = Ui, W^^Wi^ U2, W^^Wq^ Us, Wj ^Ws^ U4, 
Wq = Wio = U5, Wii = Wi2 = Uq, Wi3 = Wu = U7, W15 = Wie = Us, (3.15) 

which carried into the six constraint equations (3.6) yield six new constraints 

Us = e ui,U4 = e U2,U5 = {l + e )ui - U2, 

4K 8K /o 1 r\ 

U6 = us,U7 = e ur„us = e ui. (3.1oj 

The general expression of the rates satisfying these constraints depend on two 
independent parameters: 

w(a ge (Tn ctw crs) = (3-17) 
^ 1^1 + caw<yN + (7^ - c)(Te(Tn - (7 {^^w + 1(7n +{i - ^^'^ ■ 

Let us now consider two particular cases of this general expression. 

• The previous case where the central spin is influenced only by the North and East 
spins is a particular case of the model just considered. In order to recover the 
results (3.12) and (3.13) it suffices to impose the additional conditions 

U1 — U2 — X, Us — Ui — z, u^ — uq — z' , u-j — us — y (3.18) 

in (3.16), yielding c = in (3.17), that is (3.13). 

• If we impose the left-right spatial symmetry on the rates (3.17), i.e., if U2 — u^, 
which itself leads to u^ — u-j, then we obtain the unique solution 

/ 7^ 7 \ 

w{aaE(TN(yw(7s) = ^ U + y (^E + ^w) - i^cf{cfw + Sctat + cte) 1 , 

(3.19) 

up to the global scale a. This expression cannot be recast in an exponential form 
analogous to (3.14). 

3.4- Directed dynamics on the triangular lattice 

We now consider the two-dimensional Ising model on the triangular lattice. There 
are 128 (2^) unknown rates, corresponding to the 128 motifs where the central spin is 
surrounded by its six neighbours. Using the same method as above, we obtain the set 
of constraints between the transition rates that need to be satisfied in order for the 
stationary state measure to be Gibbsian. The result is rather lengthy and will not be 
written here. 

Let us specialize to the directed model in which the subset of influential spins are 
any three consecutive spins at 7r/3 one from the other. The sites corresponding to the 
subset of influential spins and the sites corresponding to the complementary subset of 
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neighbouring spins are thus related by spatial parity with respect to the central site. For 
this directed case, we find a unique solution, up to a global timcscalc, to the problem of 
the determination of rates leading to a Gibbsian stationary state. The solution found 
has again the exponential form 



where (7i, (72 and (73 are the spins felt by the central spin, denoted by a. 
3.5. Directed dynamics on the cubic lattice 

We finally consider the three-dimensional Ising model on the cubic lattice. Again there 
are 128 (2^) unknown rates, corresponding to the 128 motifs where the central spin is 
surrounded by its six neighbours. We find that the set of constraints between the 
transition rates that need to be satisfied in order for the stationary state measure 
to be Gibbsian are the detailed balance relationships. In other words, for the three- 
dimensional Ising model on the cubic lattice, the only possible dynamics leading to the 
Boltzman-Gibbs distribution at stationarity is the symmetric one. 

The method used in this section is unfortunately cumbersome when applied to 
lattices of increasing coordination number. The next section provides an alternate 
approach. 

4. Lattices with given coordination number 

In this section wc introduce a complementary viewpoint on the general problem of the 
existence of rates leading to a Gibbsian stationary state, when the dynamics is directed. 
We proceed in two steps. First, wc check the results of section 3 for exponential 
rates: assuming the form of the rates, we show that they satisfy, or do not satisfy, 
the master equation at stationarity, depending on the lattice considered. Secondly, 
generalizing these observations, we give a method for the determination of which lattice, 
of given coordination number, is able or unable to sustain exponential rates, for directed 
dynamics. 

4.1. Check of the results of section 3 

The idea stems from the observation that, once the form of the rates satisfying the 
master equation (2.5) have been determined by the method of section 3, it should be 
an easy matter to check the reverse, i.e., that the stationary master equation (2.5) is 
satisfied. For example, for the directed Ising chain with rates (2.28), eq. (2.5) leads to 
the following condition, for a given fixed configuration C, 



w{aaia2 ... ere) = e 



-2Ka{ai+a2+cr3) 



(3.20) 



—2K(Ji(Ji+\ _ 



(4.1) 



which indeed holds, due to the translational symmetry of the system. 
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We now consider the general case of a lattice of given coordination number z. We 
assume that the rates have the exponential form of eqs. (2.28), (3.14) and (3.20), 

«;(Ci|C) = e-^'^'^^'^S (4.2) 

where we denoted by hf the field felt by spin i due to the subset v^{i) of its influential 
neighbours. Likewise we denote by the field due to the complementary subset v~{i) 
of neighbouring spins. We only consider the case where the subsets v^{i) and v~{i) are 
related by spatial parity with respect to the central spin. With the hypothesis (4.2), 
the master equation (2.5) gives 

e-'^'^*^^ = e-'^"^'^'" , (4.3) 

i i 

or 

Y n e-2^"'"^=^ n e~2i^<^^"^ (4.4) 

i jev+{i) i jev-{i) 

that is, with r = — tanh 2K, 

Y n (l + ra,a,)=^ [J i^ + ra^a,). (4.5) 

This equation must be satisfied order by order in r. This leads to a set of geometrical 
constraints that a given lattice, of given coordination number, should satisfy, in order 
for the master equation to be satisfied. 

Before using the method in all its generality we first check that for the 2D square 
lattice, the 2D triangular lattice and the 3D cubic lattice, we recover the results found 
in section 3: the constraints are satisfied for the two first lattices, not for the last one, 
as we now show. 

We start with the 2D square lattice. We have hf = J{(Te + cfn) and = 
J{o'w + C5). The constraint equation (4.5) reads 

^ 1 + Tai{aE + (Jn) + t'^cteo-n = X) 1 + raiiaw + crs) + T^awos, (4.6) 

i i 

which should be fulfilled for any given fixed configuration C. At order r this imposes 
that the equation 

Y^i{^E + (yN) = Y^ii^W + (7s) (4.7) 
i i 

be satisfied. This is obviously the case because of translation invariance. At order 
the constraint reads 

i i 

which again is easily seen to hold. This completes the proof for the 2D square lattice. 

We proceed with the 2D triangular lattice, with z = Q. Let us denote by Ci, 62, 63 
the unit vectors spanning the lattice. Then, at order r, the constraint equation reads 
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where ei = — ei, and so on. Again, trivially it is possible to pair the terms of both sides 
by translation. At order the constraint is 

i 
i 

It is also easy to check that this equation is satisfied by translation invariance by pairing 
adequately the terms on both sides. At order r^, the constraint equation reads 

i i 

This equation is indeed satisfied for the triangular lattice because the following 
geometrical relationship between unit vectors holds: 

62 = ei + 63. (4.13) 

Thus the quadrilateral formed by the spins on the left-hand side is related by translation 
to that appearing in the right-hand side. This completes the proof for the 2D triangular 
lattice. 

Let us now consider the 3D cubic lattice, with z = Q. The constraint equations 
are the same as for the 2D triangular lattice, since these equations only depend on the 
coordination number and not on the structure of the lattice. In the present case, the 
constraints are satisfied at order r and r^, but not at order since (4.13) no longer 
holds. We thus confirm the result of section 3, namely that the 3D cubic lattice cannot 
sustain exponential rates for directed dynamics. 

4-2. Which lattice of given coordination number can sustain exponential rates for 
directed dynamics? 

As exemplified by the two last cases considered above, where z — the constraints 
coming from equation (4.5) only depend on the coordination number z of the lattice. 
Then, in all generality, for z given, the question is whether there exists a lattice fulfilling 
these constraints or not. 

Consider for example the case of a lattice of coordination z = S. The constraints 
at order r and are trivially satisfied, as above. At order the constraint imposes, 
between the four unit vectors spanning the lattice, equalities of the type (4.13). At 
order the constraint imposes the relationship: 61 + 64 = 62 + 63. This relationship can 
be implemented for the 3D centered cubic lattice, for example. However the constraints 
imposed at order cannot be fulfilled either by the 3D centered cubic lattice or by any 
other one. 

Some more work can convince the reader that any lattice with coordination number 
^ > 8 does not satisfy the constraints, hence cannot sustain a directed model with 
rates (4.2). 

This method assumes rates of the form (4.2). It is however plausible that if these 
rates do not produce stationary states described by the Gibbs measure, then no other 
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form is able to do so, because the exponential form (4.2) is the best fitted to be a 
solution of the master equation (2.5). We are therefore led to conclude that directed 
models with Gibbsian stationary states only exist in dimension one and two. 

5. Directed Ising-Glauber models 

In this section we investigate whether directed Ising models, augmented by Glauber 
dynamics, can exhibit a phase transition to a ferromagnetic state. 

5.1. The one- and two-dimensional directed Glauher-Ising models 

We begin by presenting the transition rates for directed Ising models with Glauber 
dynamics. Let hf denote, as above, the field felt by spin i due to its influential 
neighbours ?;+(?'). For example, for the directed Ising chain where each spin only sees 
the spin to its right, the field felt by spin i due to spin i + 1 is hf = Jo"j+i. In the 
two-dimensional case where the spin (jj sees the spins and a^;, the field felt by spin 
i is given by hf = J{(Jn + cte)- We define the Glauber rate for the process in which 
spin cTj is fiipped in configuration C = {cti, . . . , aj, . . . , cxiv} resulting in configuration 
= {cTi, . . . , —CTj, . . . , ctat}, as 

w{Ci\C) = ^(1 - aitanhphf). (5.1) 

We now consider the equation of motion for the magnetisation, rUi = (o"j), of the 
spin at site i, first for the directed Ising chain with periodic boundary conditions. This 
equation reads (where a dot indicates a time derivative) 

Thi — — mi-\- (tanh(i^'crj+i)) 

— — mi-\- tanh K rui+i, (5.2) 

{K = pj). The periodic boundary conditions, combined with translational invariance, 
imply that = m, independent of i, giving 

m{t) = m(0) exp[-(l - tanhis:)^]. (5.3) 

This shows that the one-dimensional system cannot support an ordered phase, since an 
initial magnetisation will decay exponentially to zero on time scale tcq = (1 — tanhi^)"^. 
Note that the stationary measure for this case is the Boltzmann distribution with 
temperature doubled, as can be seen from the alternate form of the rate (2.29). Therefore 
the vanishing of the magnetisation in the stationary state of the ID directed Ising model 
is not very surprising. However, a very similar analysis can be employed for the 2D 
directed Ising model defined above. 

Let us label the sites with a pair of integers these being the coordinates in 

the two spatial directions, and suppose the site i,j is influenced only the spins at sites 
{i + 1, j) and {i,j + 1). Then the equation of motion for ruij reads 

rhij = - ruij + (tanh[ir((7ij+i + (Jj+ij)]) 

= - mij + - tanh(2ii')(mij+i + rui+ij), (5.4) 
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where the last line follows from the fact that (7j + CTj+ij can take only three values, 
— 1, and 1. With periodic boundary conditions, translational invariance implies that 
ruij — m, independent of the site indices giving 

m{t) = m(0) exp[-(l - tanh2X)i]. (5.5) 

It follows that the magnetization again vanishes in the stationary state. 

One can also treat the case of open boundaries, where translational invariance no 
longer applies. In one dimension, for example, the right-most spin, ex at, retains it initial 
value for all t. Without loss of generality, we can fix cr at = 1. In the stationary state, 
eq. (5.2) then gives rui = (tanhi^)^"*, so the local magnetisation decays exponentially 
as a function of distance from the fixed spin ctat, and the total magnetisation per spin 
vanishes in the thermodynamic limit. 

5.2. Cayley Trees 

We now discuss the case of directed Ising models on Cayley trees, for the case of Glauber 
dynamics. We take the influence to be directed from the tips towards the root, i.e. each 
spin only sees the spins that are further from the root. 

We consider first a tree with branching ratio 2, i.e. coordination number 3. We 
assume that the initial condition is one in which all spins at a given level of the tree 
have an equivalent initial condition, e.g. all spins up, or all spins either up or down with 
probabilities that are the same for all the spins at that level. Then the Glauber equation 
for the magnetisation rUn of a site at the nth level of the tree, counting the tips as level 
1, is 

rrin = -rrin + (tanh K{S[''~^'> + 4"""^^)), (5-6) 

where, for a given spin at level n, Si"'"^^ and 82^'^^ are the two spins at the {n — l)th 
level that influence the given spin. Then, as for the d — 2 case, one has immediately 

rrin — —nT-n + tanh(2ir) m„_i. (5.7) 

In the stationary state, this gives, = tanh(2i^) mn-i. Under iteration, the site 
magnetisation is driven to zero, i.e. the interior of the tree remains unmagnetised, at 
any non-zero temperature, even when the spins at the tips are completely aligned (and 
recall that, if initially aligned, they will stay aligned forever, as nothing can change the 
tip spins in the directed model). We conclude that there is no phase transition to a 
ferromagnetic state in the directed Cayley tree with branching ratio 2. 

It is instructive to contrast this with the properties of the undirected Cayley tree. 
In this case we will partially align the tip spins by applying a magnetic field h to these 
spins, and ask whether, in the equilibrium state, the induced magnetisation propagates 
into the interior. It is trivial to trace out the tip (i.e. level 1) spins to obtain a new 
model in which the new tip spins (now level 2) experience a field h' given by 

H' = q tanh"^ [tanh H tanh K] , (5.8) 
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where q is the branching ratio, H — h/T, and K — J/T, with T the temperature, now 
a physical variable. This recurrence relation has a trivial fixed point at H — 0, which 
is stable for qtanhK < 1. In this regime, the interior of the tree is unmagnetised. 
For qtanhK > 1, however, the recurrence relation has a stable non-trivial fixed point, 
and spins in the interior of the tree have a non-zero mean value. These two regimes 
are separated by a critical coupling Kf, given by tanhi^c = 1/9- We conclude that the 
undirected model has a phase transition for all > 1, but for g = 2 the directed model 
does not. 

To explore further the directed models we consider the case q — 3. Since each spin 
at the nth level is influenced by three spins at the {n — l)th level, the analogue of eq. 
(5.6) is 

rhn = -m„ + (tanhX[>Si"-'^ + S^-''^ + 4""'^]>, (5-9) 

Exploiting the fact that {Si -|- 5*2 + Ss) can only take four values (3, 1, —1, —3), and its 
odd parity under Si — > —Si for all i — 1, 2, 3, one can write 

(tanhi^[5{"-') + 5^"-') + S^-''^]) 

= 3Amn-i - B{S[''-'^ St~'^ St''^), (5.10) 

where 

A = [tanhis: + tanh(3ii')]/4 , 

B = [3tanhX - tanh(3is:)]/3 . (5.11) 

It is clear that spins at the same level are uncorrelated, since they are influenced by 
disjoint sets of spins at lower levels. It follows that S's""^'* ^'3""^^) = In 

the stationary state, therefore, eq. (5.9) reduces to 

m„ = 3A m„_i - B ml_^ . (5.12) 

This recurrence relation has two flxed points: a trivial one m* — 0, which is stable 
for A < 1/3, and a nontrivial one m* — [{3A — 1)/BY/'^, which is stable for 
A > 1/3. There is a critical coupling Kc obtained from the condition A — 1/3, 
i.e. tanhK + tanh(3X) = 4/3, with solution = 0.4753269. We deduce that the 
magnetisation deep in the interior of the tree vanishes for K < and is nonzero for 
K > Kc, i.e. there is a phase transition at K = Kc- 

Higher values of the branching ratio q can be treated in a similar way. In each 
case one finds a phase transition at a critical coupling value Kc, which decreases with 
increasing q. For q — 4:, for example, one finds Kc = 0.3102182. 

It is interesting that the smallest integer q for which there is a phase transition on 
the undirected tree is g = 2, whereas on the directed tree it is g = 3. This is consistent 
with our observation that the 2D Ising model on a directed square lattice, in which 
each spin is influenced by only two other spins, has no phase transition and raises the 
question as to whether the 2D Glauber-Ising model on a triangular lattice, with each 
spin influenced by (say) three other spins, would have a phase transition. 
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6. Discussion and Summciry 

In this paper we have addressed two open questions. The first, is under what 
circumstances the stationary states of directed Ising models have the conventional 
Boltzmann-Gibbs measure. In such cases, there will be a conventional ferromagnetic 
phase transition identical to that of the undirected Ising model, for space dimensions 
d > 2. We have shown that there is a unique set of rates for which the Boltzmann- 
Gibbs measure describes the stationary state in dimension d = 1 and d = 2, both for 
the square and triangular lattices (where in the latter models each spin sees half of its 
neighbours). We have also shown that for the cubic lattice in dimension d = 3, the only 
set of rates satisfying the constraints of stationary Gibbs measure are those satisfying 
detailed balance. In one and two dimensions it is actually quite simple to demonstrate 
that the rates (4.2) do lead to stationary states with the Boltzmann-Gibbs measure by 
simply substituting this measure into the stationary master equation. The advantage of 
the present approach is that it shows that these rates are unique. We have generalized 
this procedure and used it as a method for the determination of which lattice is able 
to sustain transition rates for directed dynamics with stationary Gibbs measure. We 
confirmed that the 2D triangular lattice is able to do so, while the 3D cubic lattice is 
not. More generally we claim that lattices of coordination number z > 8 cannot be 
directed, if their stationary measure is Gibbsian, and conclude that directed models 
with such measure only exist in dimension one and two. 

The second, related question, concerning the existence of phase transitions in 
directed Glauber-Ising models, was inspired by the numerical studies presented in [5], 
where it was demonstrated that a directed Glauber-Ising model on a square lattice, with 
each spin seeing only its North and East neighbours, does not exhibit a phase transition 
to a ferromagnetic state. Here we have provided a simple analytic proof that there is 
no transition. We have also shown that directed Gayley trees exhibit a phase transition 
when the branching ratio q satisfies q > 3. 

A number of open questions remain. A completion of the proof of the absence of 
rates producing a Gibbsian stationary state for directed lattices of coordination z > 8 
would be desirable. It would be worth studying the 2D hexagonal lattice and see whether 
we can apply the method of sections 2 and 3 to this case. It would also be interesting to 
compare the dynamical properties of the directed and undirected models studied in the 
present work in one and two dimensions, when both share the same Gibbsian stationary 
measure. 

Secondly, the results on the Gayley tree suggest that the coordination number plays 
an important role, and raises the question as to whether a 2D directed Glauber-Ising 
model in which each spin sees more that two neighbours (e.g. a square lattice in which 
each spin sees its North, East and West neighbours, or a triangular lattice) could have a 
phase transition. Another obvious question concerns the directed 3D model, with each 
spin seeing its "North" , "East" and "Up" neighbours. Does this model, with Glauber 
dynamics, have a ferromagnetic phase transition? The data of ref. [5] suggest not. 
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Unfortunately, an analytic study for this case does not seem straightforward. 
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Appendix A. 

In this appendix we give the form of the rates for the 2D directed Ising models on the 
square lattice with Hamiltonian defined with two coupUng constants, Ji, J2- 

For the model where each spin sees its North and East neighbours, eq. (3.13) is 
replaced by 

ot 

w{(T(TE(TNcrwcrs) = ^ (1 + 7i72 o-fiO-JV - cr{jiaE + 720'Af)) , (A.l) 

where 712 = tanh2fri 2, which, after fixing the global timescale by the choice a — 
2cosh2i^i cosh2i^2, yields the unique solution: 

w{a ge ctn ctw cts) = e-2'^('^i'^^+-^2''^). (A.2) 

For the model where each spin sees its West, North and East neighbours, eq. (3.19) 
is replaced by 

CK / 7l72 

w{a Ge (Tn (yw crs) ^ 2 V ~2~ ^^^^^ ^ ~ l^ili^^w + 272(7jv + 7i(7e) 
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